Dynamical Self-Quenching of Spin Pumping into Double Quantum Dots 
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Nuclear spin polarization can be pumped into spin-blocked quantum dots by multiple Landau- 
Zener passages through singlet-triplet anticrossings. By numerical simulations of realistic systems 
with 10 7 nuclear spins during 10 5 sweeps, we uncover a mechanism of dynamical self-quenching 
which results in a fast saturation of the nuclear polarization under stationary pumping. This is 
caused by screening the random field of the nuclear spins. For moderate spin-orbit coupling, self- 
quenching persists but its patterns are modified. Our finding explains low polarization levels achieved 
experimentally and calls for developing new protocols that break the self-quenching limitations. 

PACS numbers: 73.63.Kv, 72.25.Pn, 76.70.Fz 



Double quantum dots (DQDs) are promisingplatforms 
for spintronics [l| and quantum computing [2-5]. For 
qubits encoded in singlet (S) and triplet (T) states of 
spin blockaded DQDs [6|, the hyperfine coupling of elec- 
tron spins to the nuclear spin reservoir is critical. Al- 
though electron spin relaxation caused by this coupling 
is destructive, a properly controlled nuclear polarization 
is an efficient tool for performing rotations of S-T qubits 
T is the zero component of the electron spin triplet 
(Tq,T±). A widely discussed approach for pumping nu- 
clear spin polarization into a DQD is based on multiple 
Landau-Zener (LZ) passages across the S — T+ anticross- 
ing P-[l3l| (T + is the lowest energy component of the 
triplet T in GaAs). fn the absence of spin-orbit (SO) 
coupling, angular momentum conservation requires that 
each transformation of the S state into the T + state is 
accompanied by the net transfer of one quantum unit 
of angular momentum to the nuclear subsystem. Mul- 
tiple S — > 7+ passages increase the nuclear spin polar- 
ization, but experimental data show that spin pumping 
typically saturates at a surprisingly low level of about 
1% and the origin of this puzzling behavior remains 
unknown. Higher levels of the nuclear polarization dif- 
ferences ("gradients") between the two dots were only 
achieved by using feedback loop schemes [TBI ]. 

In this Letter, we show that nuclear spin pumping pro- 
duces dynamical screening of both the SO coupling and 
the random hyperfine (Overhauser) field controlling the 
width of the S — T + anticrossing as well as the efficiency 
of spin pumping, while the detailed patterns differ. The 
screening results in one of the nuclear spin configura- 
tions with a vanishing anticrossing width, v — > 0. Hence, 
the probability of the S — > T + transition and the an- 
gular momentum transfer inevitably vanish, resulting in 
quenching of spin pumping. We call this dynamical self- 
quenching of spin pumping into double quantum dots 
borrowing the word "self-" from the theory of polarons, 
where self-trapping implies a joint evolution of the elec- 
tron and phonon subsystems jig . 

As applied to hyperfine coupled systems, this conclu- 
sion appears to agree with the concept of dark states en- 



visioned in Rcf . [17[ and further discussed in Ref . [lOj ; the 
latter paper is mostly concerned with the building of gra- 
dient fields. However, the existence of dark states has no 
direct experimental confirmation yet. On the theoretical 
side, the patterns of highly-nonlinear coupled electron- 
nuclear dynamics that might bring systems including mil- 
lions of nuclear spins into such states remain unclear. To 
resolve the problem, we performed large scale numeri- 
cal simulations for realistic systems. Our procedure (i) 
evaluates the coherent precession of the coupled electron 
spin and about 10 7 nuclear spins subject to an external 
magnetic field during a large number (up to 10 5 ) of LZ 
sweeps through the S—T + anticrossing and (ii) computes 
the building of the nuclear polarization during each LZ 
sweep. The calculations unveiled the gross features of 
the self-quenching process. Among our results, the fol- 
lowing are of special importance: (i) self-quenching sets 
in under generic conditions, (ii) spin-orbit interaction is 
dynamically screened despite the violation of the angular 
momentum conservation, (iii) durations of S — T+ pulses 
have a critical effect on the quenching dynamics, and (iv) 
dynamical screening is robust with respect to moderate 
noise levels. 

We consider two electrons in a double quantum dot 
that can be in singlet or triplet states and represent 
the orbital part of the wave function as ips( r i, r 2) or 
Vr(ri, r 2 ). The electrons arc coupled to the nuclear spins 
via the hyperfine coupling Hamiltonian 



H hf = V s ^A 



e 

j£\ m=l,2 



s(m)£(Rj-A - r m ) , (1) 



where s(m) = <x(m)/2 are the electron spin operators 
in terms of the Pauli matrices <x, m = 1,2 enumerates 
electrons, r m are electron coordinates, lj\ arc nuclear 
spins, A enumerates nuclear species and j lattice sites 
RjA; A\ are hyperfine coupling constants for the species 
A, and V s is a volume per unit cell. We consider GaAs 
that has three spin 1 = 3/2 nuclear species ^Ga, 7 ^Ga, 
and 75 As. All GaAs parameter values are known 0, [HI, 
Hil l and listed in Ref. (joj . The electrons and nuclear 
spins are subject to an external magnetic field which is 
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aligned along the z-direction. 

Before presenting our numerical results, we review 
how electronic Landau-Zener sweeps influence nuclear 
spins [13j ]. The hyperfine interaction of Eq. ([1]) cou- 
ples two electrons to a large number of nuclear spins. 
We account for the effect of this interaction by using 
a semi-classical Born-Oppenheimer aproach so that the 
slow nuclear spins produce a coupling between the sin- 
glet and triplet electronic states. In turn, the nuclear 
spins are driven by the electron dynamics controlled by 
time-dependent gate voltage variations. The variations 
causing Landau-Zener transitions occur in the interval 
— Tlz < t < Tlz and they are repeated many times af- 
ter waiting for a time T w , where T w ^> Tlz- During 
the waiting time T w , the nuclear spins are only affected 
by the external magnetic held and not by the electrons. 
During the LZ sweeps, the voltage changes also induce 
relative shifts of the electron singlet and triplet energy 
levels driving passages of the system through the S — T + 
anti-crossing, Fig. [Ha). Restricting the discussion to its 
vicinity and disregarding contributions of the (Tb,T_) 
spectrum branches, the coupled equations for singlet and 
triplet amplitudes cs{t) and CT + (t) are (h = 1) 



id i 



cs 

Ct. 



e s v+ 

V~ e T+ - V 



cs 

CT + 



(2) 



where es(t) and er + (i) are electronic energies controlled 
by the gates. The off-diagonal matrix elements of Eq. 
([2|), v ± = v^+v^, include contributions from the nuclear 
spins generated by the hyperfine coupling 



A j£\ 



(3) 



and from the spin-orbit coupling [21(. The diag- 
onal contribution 77 = rjz + ffo includes the Zeeman 
energy of the T + state in the external magnetic field 
B = Bz and the Overhauser field of the nuclear polar- 
ization, i] n = -V a J2\ A x J2je\ Cjxljx wnic h is defined 
solely by the coupling of the triplet component of the 
electron wave function to the longitudinal nuclear spin 
polarization. The singlet-triplet coupling constants are 
pjx = J dripg(r, Hj\)ipT(r, Rja) an d the coupling con- 
stants in the T + state are Q\ = J dr |V>t(i*, • 

Between the Landau-Zener cycles, during every wait- 
ing period of duration T w , the nuclear spins precess in the 
external magnetic field. It is assumed that the time-scale 
Tlz for the LZ transition (induced by rapid changes in 
the gate voltages) is much shorter than the nuclear pre- 
cession times t pr for the species in the external magnetic 
field, <7iAs, i69Ga, and tjiQa,, respectively. Because of 
the large number of nuclei in the dot, N ~ 10 6 — 10 7 , the 
changes AL^ acquired by individual spins during a single 
sweep are minor. We evaluate them by integrating the 
equations of coherent nuclear dynamics dlj/dt — Aj x L, 
over the Landau-Zener transition time —Tlz < t < Tlz- 



Here Aj are the Knight fields following from Eq. ((T|). 
Finally, AIj\ = Tj\ x where 

= -V s AxPjx (Pvy + Qv x ) /(2v 2 ), (4a) 



= VsAxPjx (Pv x - Qv y ) /{2v 2 ), 
T] x = V s AxQxR/(2v), 



(4b) 
(4c) 



with v=\v ± \ = ^{vl+vl)/2 and v ± = (v x ± iv y )/y/2. 
In these expressions, < P < 1 is the S-T + transi- 
tion probability and an unbounded real number Q is the 
shake-up parameter defined as [l3j 



P + iQ 



-i2v~ 



-T LZ 



dt c s (t)c* T+ {t) 



(5) 



in terms of the singlet (triplet) amplitudes, and R = 
2v dt I c t + (t) I accounts for the Knight shift due 

to the electron spin in the triplet T + state during the 
Landau-Zener transition. 

It follows from Eq. [2] and Eq. [5] that P is equal to 
the change of the occupation of the singlet state |cs(t)| 2 
during the sweep and therefore coincides with the LZ 
transition probability. The change AI Z in the total lon- 
gitudinal magnetization I z equals 

P Q 

AI z = -—^(v v+ + v + v n )-i—^(v v+-v+v n ). (6) 

When = 0, the second term in Eq. |6] vanishes and 
AI Z = — P as required by the conservation of the angu- 
lar momentum. In the same scenario, Q ^ and de- 
scribes the angular momentum transfer inside the DQD 
due to the shakeup processes induced by the LZ pulses 



13| . When 7^ 0, Q also mediates the angular mo- 



mentum leakage from the DQD due to the spin-orbit 
coupling. Because the integral for P converges fast at 
the scale of t ~ l/l^l while the integral for Q diverges 
as In Tlz for linear LZ sweeps, Q is typically large, espe- 
cially for P rj 1, and deeply influences the self-quenching 
process. 

Using the amplitudes (cs(t),CT + (t)) found from solv- 
ing Eq. ^ in combination with the dynamical equations 
for nuclear spins of Eq. Q makes our approach com- 
pletely self-consistent. During a single LZ sweep, es(t) 
and €r + (t) change fast while staying close to the anti- 
crossing, and rj and remain practically constant. The 
singlet wave function of the DQD in the S — T+ anticross- 
ing point can be expressed as ips = cos v 1P02 + sin v ipn, 
where ipQ2 and ipn are the singlet wave functions with 
both electrons on the right dot and two electrons equally 
distributed between the dots, respectively. The mixing 
angle v is controlled by B and detuning and was chosen 
as v = 7r/4. 

Our simulations included about 10 7 nuclei up to those 
with a hyperfine interaction strength of only 1% of the 
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maximum one, which allowed us to account for the elec- 
tronic density heterogeneity and the spin polarization 
transfer from the interior to the periphery. Our algo- 
rithm allowed performing calculations for different Tlz 
with the identical initial distribution. 

For each sweep, the DQD is first set in its eigcnstate 
at t = —Tlz that is close, but not identical, to the sin- 
glet (0,2) state with both electrons localized at the right 
dot. Then a change in the gate voltages drives a (partial) 
transition to the triplet (1,1) state with electrons shared 
between both dots. Finally the electronic system is re- 
set in its initial state. We assume that the LZ transition 
time Tlz are much shorter than the nuclear precession 
times t pr and compute the change in the direction of each 
of the nuclear spins during every sweep numerically, as 
described by Eqs. (HI-©. Between consecutive sweeps, 
repeated with a period of the waiting time T w , electrons 
are in the singlet state and do not interact with the nu- 
clear spins that coherently precess in the external field. 
We choose realistic parameters for a parabolic DQD of 
a height w = 3 nm, size t = 50 nm, and interdot sep- 
aration d = 100 nm; magnetic field Ti=10 mT. All the 
results presented below were found for the same initial 
configuration of nuclear spins, but we have checked that 
they are representative for generic initial configurations. 

For the shape of the S — > T + pulses, we used the 
LZ model with e s (t) = e max t/2Ti,i and — 77 = 
-(emaxV^Lz) - (»7 - »7i), where 77, is the initial polar- 
ization r\i = rj(t = —Tlz) and e max = 2.5 meV, which is 
larger than the typical S-T + coupling. To avoid trivial 
quenching due to the shift in 77 caused by the accumulat- 
ing polarization, the electronic energies were renormal- 
ized after every 100 sweeps to keep 77 — rji « 0. As a re- 
sult, the center of the sweep was permanently kept close 
to the anticrossing point. Such a regime can be achieved 
experimentally by applying appropriate feedback loops. 
SO coupling in DQDs is device specific, and in GaAs it 
changes from weak to moderate, so we consider it both in 
the limit of no SO coupling [71 1 2 21 and with SO coupling 
of a reasonable magnitude [23j, |24j . 

Fig. QJb) plots the change in the total nuclear polar- 
ization, A/ z , as a function of the number of sweeps n, 
for v so = and four transition times Tlz; the differ- 
ence in the parameter values of all three nuclear species 
is taken into account. The evolution of AI Z typically 
saturates within 3 x 10 4 sweeps. The saturation proves 
the self-quenching of the transverse nuclear polarizations 
that controls the singlet-triplet coupling shown in Fig. 
E^c). It vanishes after a number of LZ transitions and this 
typically happens faster for longer LZ durations Tlz (a 
larger transition probability) , but more complicated pat- 
terns of subsequent revivals of the i;„ can also be seen. 
Volatile dynamics of seen in Fig. Qlc), with multiple 
maxima and minima, is typical of multi-specie systems 
because of the different spin precession rates of differ- 
ent species. However, finally the nuclear subsystem self- 



synchronizes in one of the states in which it decouples 
from the electron spin qubit, and this is our first central 
result. The contribution of each of the species to van- 
ishes identically, at each instant; v„ = persists even 
after the LZ pumping is interrupted (not sho wn) . This 
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result resembles the "dark states" of Refs. 

With N ~ 10 6 nuclei in the DQD, the initial fluctu- 



ations producing is 



N 1 / 2 - 10 3 . For P ~ 1, one 
expects that at least n ~ 10 3 pulses are needed for bal- 
ancing it. The typical number of pulses to establish self- 
quenching of about n ~ 10 4 of Fig. [T^b) is one order of 
magnitude larger, which can be attributed to the high 
volatility of the process and the fact that the LZ proba- 
bility P in each cycle is less than one. 
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FIG. 1: Nuclear dynamics in absence of SO coupling as a 
function of sweep number n; difference in the parameters of 
the three nuclear species is taken into account. Resonant 
waiting time T w = t75As = 13.7/xs; qualitative pattern do not 
depend on this specific choice, (a) Landau-Zener passage of a 
singlet S through a S — T+ anticrossing. Energy levels (black) 
and evolution of S into entangled S and T+ states (red), (b) 
Change in the nuclear spin polarization AT and (c) hyperfine- 
induced singlet-triplet coupling |iv |. Color codes: Tlz = 10 
ns (red), 20 ns (green), 40 ns (blue), and 80 ns (orange). 



Next, we consider the effect of the SO coupling. To 
illustrate the main qualitative result, we ascribe to all 
nuclei identical parameters found by averaging over the 
three GaAs species [2(j and consider strictly resonant 
pumping, T w = rj GaAs = 10.7/xs, t Ga As being the nu- 
clear precession time 25| . We use realistic values of the 
0, 31, and 62 neV [26| (and choose 
The LZ transition probability 



SO coupling, v s 
it to be a real number), 
shown in Fig. [21a) vanishes at large n, hence, the pump- 
ing is self-quenched also in this case. At this point, it 
is crucial to note that all our numerical results for 
and AI Z are plotted at multiples of the waiting time 
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T w . Hence, sclf-qucnching in a resonant (T w = tGaAs) 
SO coupled system is achieved through — > —v so 



■) at every multiple of the Larmor pe- 



riod, and nuclear polarization screens the SO coupling, 
Fig. [2] This screening of SO coupling is out next cen- 
tral result. However, in contrast to the case of no SO 
coupling, between the sweeps (not shown), the matrix 
elements v^(t) change harmonically with the amplitude 
v so and a period toaAs, v n -> v so cos (2-Kt / to&As) , where 
the time t = nT w = nfoaAs at the LZ sweep number n. 
Not surprisingly, while | reaches its TLz-mdcpcndent 
limit, the change in the polarization, AI Z , depends on 
Tlz. Fig. [^a) shows that P is large, P ~ 1, and fluc- 
tuates fast with n. The mechanism of fast dynamics is 
unveiled by the high magnitude of the shakeup parame- 
ter Q ~ 20 ^> 1. While P describes pure injection of the 
angular momentum, Q describes its redistribution due to 
shakeup processes and the SO coupling [l3j ■ Remarkably, 
it is seen from Eq. ([6]) that the effect of the Q term on 
AI Z vanishes in two important limits, when v so = and 
= — v so - Fig. 3(c) demonstrates that self-quenching 
sets in sharply for large Tlz, and the fluctuational phase 
lasts longer for larger v so values; always saturate at 
— v so . For SO coupled systems, the change AI Z in the 
magnetization is non-monotonic in n and shows no reg- 
ular dependence on v so , Fig. OJd). 
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FIG. 2: 

Nuclear dynamics for variable spin-orbit coupling v so driven 
by long pulses Tlz = 160 ns; single-specie model with 
Tw = tGaAs- (a) Landau-Zener probability P, (b) shakeup 
parameter Q, (c) hyperfine-induced singlet-triplet coupling 
|«n |, and (d) change AI Z in the total nuclear polarization; in 
(c), dashed lines mark spin-orbit coupling. Color codes: 
v so = (red), 31 neV (blue), and 62 neV (green). 

To test how robust our results are, we modeled the in- 
fluence of noise by adding a random magnetic held along 
the z-direction for each nuclear spin so that the nuclear 
spins acquire an additional phase of 2irrj\T w /T during 
each waiting time between LZ sweeps, where rj\ are ran- 
dom numbers in the interval from -1 to 1; r is the noise 
correlation time (see Ref. (2^) for more details). 



Fig. [3] demonstrates the effect of the noise that ran- 
domizes the phases of the precessing nuclear spins. It 
displays the magnitude of the hyperfine matrix element 
at every multiple of the waiting time T w for no noise 
and for two noise levels. In all cases, approaches the 
value sa —v so at n < 50000. While in the absence of the 
noise the saturation of is exact at all multiples of the 
waiting time T w , = — v so (blue curve), noisy systems 
experience slight fluctuations near this value, see espe- 
cially the red line which only saturates at n « 5 x 10 4 . 
We conclude that for moderate noise levels, the dynami- 
cal screening is robust, and this is our final central result. 
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FIG. 3: Effect of noise on nuclear dynamics. The LZ sweep 
duration of Tlz = 80ns; single-specie model with T w = fcaAs- 
Hyperfine-induced singlet-triplet coupling \v^\ for spin-orbit 
coupling Vso = 62 neV and increasing levels of transverse 
noise. Black dashed lines show SO coupling. Color codes: 
r/tGaAs = oo (blue), 5000 (green), and 2500 (red); tGaAs is 
the nuclear spin precession time and r is the noise correlation 
time. 



The above results prove that dynamical sclf-qucnching 
is a rather generic property of the GaAs type DQDs 
pumped b y m ultiple passages through the S — T + an- 
ticrossing [27[. In a quenched state, the electron spin 
qubit becomes screened from the randomness of the nu- 
clear spin bath, and therefore its decoherence by nuclei 



30J is expected to be suppressed; trapping the qubit 
into the quenched state can be checked by spin splitter 
technique 3l|, |32j . A quenched qubit can be operated by 



short pulses applied to the gates due to the different de- 
pendence of and v so on the shape of the electron wave 
functions. This subject requires further consideration. 

In conclusion, self-quenched states produced by sta- 
tionary pumping of a spin blockaded double quantum 
dot by multiple passages through the S — T + anticross- 
ing decouples the electronic qubit from the nuclear spin 
bath. In such states the dynamical nuclear polarization 
screens both the initial random Overhauser field and a 
moderate spin orbit coupling typical of GaAs quantum 
dots. Self-quenching unveils the origin of the low spin 
pumping efficiency encountered in experimental studies. 
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EPAPS SUPPLEMENTARY MATERIAL FOR "DYNAMICAL SELF-QUENCHING OF SPIN PUMPING 

INTO DOUBLE QUANTUM DOTS" 

In this supplementary, we provide details of the parameters used in the numerical calculations, describe the effective 
single-specie model, show how nuclear spins precess in an external magnetic held between Landau-Zencr sweeps, and 
discuss the noise model. 

Parameters 

GaAs has 8 nuclear spins per cubic unit cell so that the effective volume per site is V s — a 3 /8, where the lattice 
constant is a = 5.65A. When all nuclear spins are fully polarized, the Overhauser field seen by the electrons is 5.3 
T. The electron g-factor is go&As = —0.44. The other parameters reflecting the abundance, g-factors, and hyperfine 
coupling constants are listed in Table fl] 





69Ga 


71Ga 


75As 


p 


30% 


20% 


50% 


g 


1.3 


1.7 


0.96 


A (fieV) 


77 


99 


94 


I 


3/2 


3/2 


3/2 



TABLE I: Nuclear abundances p, nuclear p-factors, hyperfine coupling constants A, and nuclear spin in GaAs. 



Effective Single-Specie Model 

To illustrate the main qualitative effects of SO coupling, we also make use of an effective single-specie model. In 
the effective single-specie model, the mean nuclear ^-factor is <?GaAs = ^2\P\9\ = 1-21, the mean hyperfine coupling 
constant is AcaAs = ^2\P\A\ = 90//eV, and the mean nuclear spin is IcaAs = ~^2\P\I\ = 3/2. 

Nuclear Precession in the External Magnetic Field 

Landau-Zcner sweeps of a duration Tlz are repeated with with wait periods of T w long compared with Tlz, hence, 
7lz T w . Between consecutive Landau-Zcner sweeps, electrons are in the singlet state and do not interact with the 
nuclear spins. However, during the waiting interval T w , nuclear spins precess in an external magnetic field B applied 
along the z-direction. Then, the coherent evolution of the nuclear spins are 

Ijx (t + T w ) = i? A (t) + I* x (t) cos 4> jX ~ i* (t) sin (f>j\ , (7a) 
J* (t + T w ) = 7J A (t) + 7J A (t) cos + i? A (t) sin fa , (7b) 
I!x(t + T w ) = I jX (t), (7c) 

where t is the initial time of the interval, the transverse phase change is <pj\ = — 2-7rT w /i,\ in terms of the spin 
precession time t\ = 2irh/ gx^iB, gx = px/I\ is the p-factor for a nuclear specie A, /i> its magnetic moment and 
Hi = 3.15 x 1CP 8 eV/T is the nuclear magneton. 

Noise 

We neglected in our calculations the effects of the magnetic dipolc-dipole interaction between the nuclear spins that 
is slow at our time scale and also the possible electrical noise on the gates . The detailed study of these effects is 
beyond the scope of our paper. 

We use a simple model to consider the influence of such noise by phcnomcnologically adding a random magnetic 
field along the z-direction for each individual nuclear spin so that the accumulated phase in Eq. [7] changes, <pjx — > 
'f'jx = + <t>jx lsc = ~27rT w (l/iA — rjx/r), where rjx are random numbers in the interval from —1 to 1 for each 
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nuclear spin. This procedure simulates a randomization of the transverse components of the nuclear spins after a 
time of the order r, and r is termed the noise correlation time. 

Quantitively, the noise averages of the phase-dependent factors in Eq. [7] become 



1 f 1 , ,eff , sin(27rT w /T 

2y_/ rC0S ^ =C ° S ^ 2*T w /r 
1 f 1 , . , off . , sin (2^T w /r) 



dr cos 0$ = cos ^ ,1 , (8) 

ieff 

dr sin (j)f x = sm ^ A ^ ^ , (9) 



so that the phase-coherent information contained in the transverse nuclear spins components are lost after a time of 
the order of the noise correlation time r. By introducing noise, any possible build-up of the matrix elements is 
therefore lost after a time of the order of the noise correlati on time r. 



